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We introduce Random Matrix Models for the Hermitian Wilson-Dirac operator of QCD-like the- 
ories. We show that they are equivalent to the e-limit of the chiral Lagrangian for Wilson chiral 
perturbation theory. Results are obtained for two-color QCD with quarks in the fundamental rep- 
resentation of the color group as well as any-color QCD with quarks in the adjoint representation. 
For Nc = 2 we also have obtained the lattice spacing dependence of the quenched average spectral 
density for a fixed value of the index of the Dirac operator. Comparisons with direct numerical 
simulations of the random matrix ensemble are shown. 
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1. Introduction 

Chiral Random Matrix Theories [1,2] have been successful in describing lattice QCD Diiac 
spectra on the scale of the eigenvalue spacing. It was shown that they are equivalent to the £-limit 
of QCD which is given by the £-limit of chiral perturbation theory [3]. Recently, Random Matrix 
Theory was extended to include discretization effects of both the Wilson [4] and the staggered 
Dirac operator [5]. They are equivalent to the £-limit of Wilson chiral perturbation theory [6, 7] 
and staggered chiral perturbation theory, respectively [8]. 

Starting from the chiral Lagrangian of Wilson chiral perturbation theory in the microscopic 
domain, exact results were obtained for the spectral density of the Hermitian Dirac operator both 
for the quenched case [4] and the case of dynamical quarks [9]. The spectral density of the non- 
Hermitian Wilson Dirac operator, could only be accessed by means of powerful random matrix 
techniques. The results for dynamical quarks show that depending on the value of the low-energy 
constants either an Aoki phase or a first order scenario is possible [10]. Lattice results for the 
eigenvalue density [11, 12, 13] have been compared successfully to the exact results [14, 9, 15, 10, 
16] for the spectral density in the microscopic Umit. 

Recently, a great deal of attention has been focused on the conformal limit of QCD and QCD- 
like theories. Both the two-color theory and the any color adjoint theory are relevant for technicolor 
theories [17]. The advantage of SU(2) theories is that they require less fermions for achieving 
conformality than SU(3) theories. Furthermore, the SU(2) theory with two adjoint fermions is 
relevant for minimal walking technicolor theories [18]. Studies of this theory have been performed 
for unimproved Wilson fermions [19] and for the analysis of the conformal window, it would be 
useful to have a better understanding of the discretization errors. 

As is the case in the continuum theory, also at non-zero lattice spacing, there is a one to one 
correspondence between patterns of chiral symmetry breaking and the anti-unitary symmetries of 
the Dirac operator [2]. We therefore can distinguish three distinct classes. QCD in the fundamental 
representation with three or more colors is the case without anti-unitary symmetry. When we 
have an anti-unitary symmetry, [T, D] =0 for the Dirac operator D, then there are two different 
possibilities. Either = 1 or = — 1. In the first case it is always possible to find a gauge 
field independent basis for which the Dirac operator is real. This is the case for QCD with two 
colors in the fundamental representation where the chiral symmetry breaking pattern is SU(2Nf) — >■ 
USp(2Nf). In the second case it is possible to find a gauge field independent basis in which the 
matrix elements are expressed as self-dual quaternions. This is the situation for QCD in the adjoint 
representation where the pattern of chiral symmetry breaking is SU(2Nf) S0(2Nf). 

The goal of this paper is to study the effect of a finite lattice spacing on the low lying Dirac 
eigenvalues and to understand the behavior of the spectral gap of D5 + nijs = 75 (Dw + m) at finite 
quark mass as a function of the lattice spacing. In mean field theory, closure of the spectral gap 
will serve as an order pai^ameter for the onset of the Aoki phase [4]. 

This paper is organized as follows. In section 2 we introduce a Wilson Random Matrix Model 
for SU(2) with fundamental quai^ks and QCD with adjoint fermions. In particular, we consider the 
Nf flavor partition function as well as the partially quenched partition function. In section 3, we 
compare analytical results with Monte Cai^lo data of the Random Matrix Model. 
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2. Random Matrix Theory 

The random matrix for the Hermitian Wilson Dirac operator proposed in Refs. [4, 20] is given 



by 



D5= \ t h (2-1) 



where A and B are Hermitian nxn and (« + v) x (« + v) matrices, and W is a complex nx (n + v) 
matrix. All matrix elements are distributed according to a Gaussian probability distribution. The 
Random Matrix Theory for j8 = 1 and jS = 4 is obtained by simply choosing real or quaternion 
matrix elements. 

For a = 0, the Wilson-Dirac operator has v generic zero modes in accordance with the Atiyah- 
Singer index theorem. At finite a, one can define the index of the Dirac operator for a fixed gauge 
field configuration through spectral flow lines or equivalently by [21] 

v= £ sign{{k\r5\k)), (2.2) 

where the above sum is restricted to the real modes since the eigenfunctions corresponding to 
complex modes have zero chirality. 

The partition function of D5 with A^f flavors is given by 



yRMT.v 



J dD5det^f{D5 + mY5+z)P{D5), (2.3) 



where ^(Dj) is the probability distribution of the matrix elements of D5. In the microscopic limit 
where the combinations m = 2mn , z = 2zn and cf = are kept fixed as n — )• 00, the Random 

Matrix Theory reduces to the £-limit of Wilson chiral perturbation theory, 



Z]^^ = J dn{U)det''Uexp 



(2.4) 



with U e U(2A^/) /Sp(2A^/) and k: = v/2 for jS = 1 while U G U(2A^/) /0(2Nf) and k: = V for j3 = 4 
as is the case for t? = [22] . 

The full partition function at fixed vacuum angle 6 is given by a sum over the Fourier compo- 
nents Zy, 

00 

Z(0) = X ^"'2v. (2.5) 

V — — 00 

In order to access the spectral properties of the Dirac operator we employ the supersymmetric 
method of RMT. The generating function for the resolvent of D5 is the partially quenched partition 
function obtained by adding an additional fermionic and bosonic quark to the Nf flavor partition 
function 



= (det^/(r3(Dw + .n)) f;^f" + ;^+^M • (2.6) 
^ ' \ det(75(Dw + m)+z')/^ 

The resolvent of is given by 



(z, m; a) = lim d.^^], = { tr ) . (2.7) 
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Figure 1: The analytical results (solid curves) 
compared to the results of the Monte Carlo 
simulation (histograms). 
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Figure 2: The distribution of the first positive eigenvalue, with 
m = , V = 0. The average position of this eigenvalue (denoted 
by the vertical bar) shifts away from the origin for increasing a. 



Its discontinuity across the real axis gives the spectral density 

py(X^,m;a) = — Im[G^(z = X^,m;a)]. 
n 



(2.8) 



In the microscopic limit the generating function reduces to a supersymmetric extension of 
of the partition function (2.4). However the integrals over the non-compact part of U are only 
convergent for imaginary a. To obtain an analytical continuation to real a we have to rotate U -^iU . 
This results in the partition function 



'Wf+l|l 



-StrM{U-U' 



X exp 



(2.9) 



The integration manifold ?7 G U(2Nf + 2|2)/UOSp(2Nf + 2|2) is the same asfora = 0. The mass 
matrix is given by M = diag(m, m,m,m) whereas the axial mass is given by Z = diag(z,z,z',z'). 
Note that for z' — z we recover the A'f flavor partition function. 

To evaluate this integral we need an explicit parameterization of U. Since it involves only four 
Grassmann variables, it can easily be evaluated by brute force. The results are quite lengthy and 
will be published elsewhere [23]. 



3. Analytical and Numerical results 

To illustrate our analytical results we compare in Fig. 1 the results obtained from (2.9) for 
V = and m = with numerical ones obtained by calculating the eigenvalues of an ensemble of 
random matrices. Surprisingly, the spectral density at zero decreases by a factor \/2 for any nonzero 
value of a. The reason for the nonuniformity of the a — )• limit is that for a 7^ the convergence 
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Figure 3: Spectral density at v = 2 and m = 0. 
Note the presence of zero modes for a = and the 
widening of the peak as we increase a. 



Figure 4: Spectral density at v = 2 and 
in ^3. For a~0 the spectrum has a gap of 
width 2m which closes for increasing a. 



of the integral is achieved through the U^-term while for a = the convergence comes from the 
U -term. The effect of diagonal blocks that comprise the Wilson term in the Random Matrix Model 
can thus be seen for arbitrarily small values of a. 

In Fig. 2 we study the distribution of the first positive eigenvalue. Apparently, the diagonal 
blocks of the Wilson Dirac operator lead to a weak repulsion of the two eigenvalues closest to zero, 
but for larger values of a, there is no repulsion away from zero. The average position of the first 
eigenvalue increases as is shown by the corresponding vertical bar perpendicular to the real axis. 

For V 7^ 0, the distribution of the zero modes is a Dirac delta function for c? = 0. As we increase 
the value of a their distribution gets broadened with a width proportional to a. This is similar to 
what happens in QCD with three or more colors in the fundamental representation. As can be seen 
from Figs. 3 and 4, at about a = 0.5 the peak due to the would be zero modes has disappeared 
almost completely. 

For a = the spectrum of D5 has a gap [—m,m], but at finite lattice spacing eigenvalues of 
tail states penetrate the gap [4]. Our results provide an explicit analytical handle on these states 
and allow us to identify the point where eigenvalues approach the center of the spectral gap and 
inversion of the Dirac operator becomes very difficult. 

For a = 0, the spectral density of the two-color theory develops a square root type of singularity 
at the edge of the gap, P5(/L^) ~ 1 /-y/ {X^y -m^ + v5(A^ -m) (see Fig. 4). For j8 = 2, on the 
other hand, the spectral density approaches a finite limit at = m . 

In Figs. 5 and 6 we show scatter plots of the eigenvalues of Dw- For a = the Wilson Dirac 
operator Dw is anti-Hermitian and the eigenvalues lie on the imaginary axis (see Fig. 5). In contrast 
to this behavior, Dw is non-Hermitian at finite a. Hence, it has a complex spectrum. Because Dw 
is still 75-Hermitian its complex eigenvalues occur in complex conjugate pairs, and at least |v| 
eigenvalues are real (see Fig. 6). The additional real modes always appear in an even number. In 
Fig. 6 we show a spectrum of Dw for v = 5 with seven real modes. 
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Figure 5: Spectram of a randomly generated matrix 
Dw with V = 5 and m = at vanishing lattice spacing, 
i.e. a = 0. 



2.5r 

2 
1.5 
1 

0.5 
ltn(z)o_ 

-0.5 
-1 

-1.5 
-2 

-2.5- 



-0.3 



; ++ + 
+ + * 
t * 



* **/y 



+ + 
t +.*+ 



-0.2 



-0.1 





Re(z) 



0.1 



0.2 



0.3 



Figure 6: Spectrum of a randomly generated ma- 
trix Dw with V = 5 and m = for a finite lattice 
spacing (a = 1). 



4. Conclusions 



We have introduced Random Matrix Theories for the Wilson Dirac operator of QCD like 
theories and have obtained explicit analytical results for the spectral density of two-color QCD. 
The analytical results for v = have been compared to Monte Carlo simulations of the proposed 
random matrix ensemble. Furthermore, numerical results for the spectral density for non-zero 
quark mass and index v have been presented. We stress that although the increase of computational 
power has allowed for lattice simulations in the deep chiral regime it is not possible to invert the 
Wilson Dirac operator when eigenvalues are sufficiently close to zero. Our results identify the 
parameter domain where such eigenvalues appear and can be potentially useful for identifying the 
parameter domain for simulations with dynamical quarks. Actually, the probability to obtain small 
eigenvalues is higher for the two-color theory than for QCD with more colors because of the lack 
of repulsion from the origin. Our analytical results can be extended to the case of arbitrary v and 
also to the case of adjoint QCD. These will be presented in forthcoming publications. 
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